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ABSTRACT 

A method for calculating a component's design 
survivability by incorporating finite element analysis 
and probabilistic material properties was developed. 
The method evaluates design parameters through direct 
comparisons of component survivability expressed in 
terms of Weibull parameters. The analysis was applied 
to a rotating disk with mounting bolt holes. The 
highest probability of failure occurred at, or near, 
the maximum shear stress region of the boit holes. 
Distribution of material failure as a function of 
Weibull slope affects the probability of survival. 
Where Weioulx parameters are unknown for a rotating 
disk, may be permissible to assume Weibull para- 

*» well as the stress-life exponent, in order 
to determine the effect of disk speed on the proba- 
bility of survival. 

NOMENCLATURE 

c material stress-life exponent 

e Weibull slope of material modulus 

L stress cycles to failure 

L i stress cycles to failure of i th component 

L o Stress cycles where 37.7 percent of 

specimens survive 

L , stress cycles to failure of total system 

L u stress cycles where all specimens survive 

s statistical fraction of specimens that survive 

a given number of stress cycles 


Fellow, ASME. 


s t total system survivability 

V stressed volume 

T « cress level applicable for failure theory 

Subscripts : 

a, b states a and b, or bodies a and b 
INTRODUCTION 

A key component in any structural design process 
is establishing a minimum acceptable service life and 
then determining the potential of component or struc- 
tural failure before that minimum life is achieved. 
Although components may be tested to failure tc deter- 
mine their functional life, this method is certainly 
impractical in terms of time requirements, the diffi- 
culty in simulating all possible loading conditions, 
and the expense involved in full-scale testing. Also, 
early-stage designs generally undergo quite a few 
iterations until a final set of specifications is 
formalized . 

Many applications, especially in the aerospace 
industry, include size and weight constraints. This 
precludes utilizing life-enhancing design features, 
such as ultraconservative safety factors, or using 
stronger, but heavier, materials. In order to realize 
the lightest, meat compact design, components are 
purposely designed for finite life (i.e., expected 
failure), with very high probability of survival 
throughout the entire duty cycle. The most efficient 
system operation then involves regularly scheduled 
component replacement, so that the predetermined time 
between replacements .s only slightly shorter than the 
anticipated life of the component. 

Often, emp ir ica 1 -ba eed design techniques are 
employed for fatigue* tolerant designs. Although such 
procedures may result in successful designs with 
respect to fatigue considerations, the designs will 
probably not be efficient, with respect to size or 


) 


cost. Also, ths urs of new mstsrisls without sn expe- 
risncs data bass will poss a problem in making dssign 
decisions. For these reasons analytical-based life 
design methodologies are required to predict component 
life and survivability without the need for extensive 
testing or field experience. 

First-order life prediction methods are generally 
based simply on yield and fatigue-limiting stresses. 
Often, they do not account for geometry and size 
effects, surface condition, duty cycles, and environ- 
mental factors. The material life data are also fre- 
quently misinterpreted. Stress-life curves are often 
assumed to be completely deterministic when in fact 
the stress-life data are presented aa statistically 
averaged values. Finally, most l L fe analyses are pri- 
marily concerned with a few critical areas in the 
component (i.e., the areas of highest stress). How- 
ever, for designs with low stress gradients areas of 
probable failure are not intuitively obvious. 

Finite-life component design requires a probabi- 
listic approach that couples operating life with an 
expected rate of survival. Weibull (1939) demonstra- 
ted a statistical analysis that was particularly 
effective in describing experimental fatigue data. 
Lundberg and Palmgren (1947) applied Weibull analy- 
sis to contact-stress, high-oycle-f at igue bearing 
problems. Grisaffe (1965) demonstrated the applica- 
bility of a Weibull-based analysis to other types of 
durability problems. Ioannides and Harris (1985) and 
Zaretsky (1967) proposed a generalized Weibul 1-based 
methodology for structural fatigue life prediction 
based on a di screte-st ressed-volume approach. 

Zaretsky et al. (1989) coupled this methodology to a 
stress field determined by finite element methods to 
predict the life and reliability of a generic, rotat- 
ing disk. They also demonstrated the applicability of 
the methodology to the des'gn process through para- 
metric studies that showed component life sensitivi- 
ties to such design variables as disk diameter and 
thickness and bolt hole size, number, and location. 
Nemeth et al. (1990) developed a computer program for 
quantifying reliability that is based on inherent 
flaws found in ceramics. T.ie program calculates the 
fast fracture reliability of microscopically isotropic 
cer .iHiC components. This method is also based on the 
component's entire stress state, not just the maximum 
stress point. 

Although these methods allow the determination of 
a total component life, it is often necessary to exam- 
ine critical areas of the cumponent in order to pre- 
dict the probability of local failure. In this manner 
designs could potentially be optimized away from known 
cn : icai areas through appropriate sizing. Therefore, 
the objectives of this investigation are (1) to extend 
the method of Zaretsky (1987) to allow for calculating 
the local probability of failure within any compo- 
nent's stressed volume as well as within the entire- 
component and (2) to demonstrate the technique on a 
generic disk by examining the sensitivity of stressed- 
volume survivability to uncertainties in the material 
propert ies . 


ANALYTICAL METHODOLOGY 

Statistical Failure Theory 

Expet imenta 1 fatigue data can often be plotted as 
a straight line on Weibull probability paper. For a 
constant stress level the number of stress cycles to 
failure (i.e., life) is plotted on the abeciesa. The 
percentile of specimens that survive at a given life 
is plotted on the ordinate. The transformation of 


stress cycles and survivability into Weibull coordi- 
nates is given by the equation 

In In - - e In (L-L ) - e In L (1) 

S u 

where S is statistical fraction of specimens that 
survive a given number of stress cycles, e is the 
Weibull slope, L is stress cycles to failure, L is 
stress cycles where all specimens survive, and L 
is stress cycles where 37.7 percent of specimens sur- 
vive. This is referred to as a three-parameter 
Weibull equation. A two-parameter Weibull equation ib 
obtained by setting * 0. The Weibull slope pro- 
vides an indication of the scatter of the statistical 
properties, with e * 1, 2, or 3.57 being representa- 
tive of exponential, Rayleigh, and Guassian distribu- 
tions, respectively. 

The number of stress cycles to failure at a given 
stress level can be related to stress cycles to fail- 
ure at new stress level by the equarion 

. . hr 


where L* is the known life at stress level 
c is the stress-life exponent obtained from 
testing. Weibull (1939, 1951) expressed the 
ity of survival as 

In i ~ T C L*V (3) 

S 

where V is the stressed component volume. By uSir.g 
Eq. (3), the siz^ effects on survivability for equally 
stressed volumes can be expressed as 

* t * « • , 

S r* S . 

After component life and survivability have beer, 
adjusted on the basis of stress levels and stres&ec 
volume, all other combinations of life anc survivabil- 
ity can be calculated. This allows valid comparison 
of relative component life at equal survivability or 
relative component survivability at equal life. 

Using Eqs. (1) to (3), Zaretsky (1987) devei.pea 
an expression for predicting life at a given surviva- 
bility of a stressed component b relative to the 
known life of a stressed component a: 


T , and 

A 

coupon 
r vbabil - 
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This equation gives the component's . i f e relative t _ 
another similar component when both have the same 
survivability. The system life can then be calculated 
from the individual component lives as follows: 



where L ? is the total system fe anJ l is tbe i” 
component life (Zaretsky, 1987). Zaretsky*' s method 
can be extended further to examine relative component 
survivability at equal life. This is given by 


S, * S 
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Equation (7) can be used to identify critical fatigue 
components within a system and to optimize noncritical 
onaa. Bacauaa a minimum lifa ia raquirad of all com- 
ponanta in tha syatam, individual componant surviva- 
bility at that life can ba aaaaaaad ralativa to a 
pradetarminad critical componant. In thia manner, 
components whose survivability is judged to ba too low 
can ba radasignad for graatar survivability. Also, 
components whose survivability is much greater than 
that of tha critical componant can ba radasignad for 
smaller size or lower weight while still maintaining 
adequate protection against failure. 

Finally, the component's survi .ability is equal 
to tha product of all the elements' survivabilities: 

s B « n s 1 < 8 > 

I ncorporating Finite Element Analysi s 

The methodology just described can easily be 
incorporated into the finite element method to deter- 
mine the life and survivability of a structural compo- 
nent. The entire component is analogous to a system. 
The finite element model of the component discretizes 
its geometry into elemental volumes. These elements 
are considered to be the base members of the "system.” 
The time to crack initiation for each element is cal- 
culated from the elemental stress levels and volume as 
descr .bed previously to determine the component's 
incipient failure time. Elemental survivability is 
also calculated to find areas that have either too low 
survivability or are overdesigned. 

A relative comparison approach is used in this 
methodology. A critical element in the* model is 
selected on the basis of a maximum stress state. The 
selected stress failure criterion can be based on the 
material used in the component. If no fatigue data 
are available, the stressed element can be arbitrarily 
assigned a life and a survivability that is used to 
normalize the life and survivability of the other ele- 
ments. This approach allows easy, qualitative compar- 
isons between designs. Also, only one set of coupon 
fatigue tests at operating temperatures is necessary 
to fix the Weibull parameters and to establish 
quantitative component lives if the stress-life 
exponent of the material is already known. However, 
if the stress-life exponent is not known, at least two 
additional sets of coupon fatigue tests will be 
necessary . 

A methodology for predicting elemental surviva- 
bility at a fixed life is illustrated in Fig. 1. The 
critical element's life and survivability is set at 
unity and 90 percent, respectively (point 1). This is 
the analysis point to which all the other elements' 
lives and survivability will be referenced. Other 
designs can also be referenced to this point by using 
the same critical volume and stress. The expected 
survivability at any number of stress cycles can be 
determined by using Eq . (1). This is shown in the 

figure by drawing a line through the reference point 
with a slope equal to e, the material’s Weibull 
slope . 

Next, an element's survivability is adjusted 
based strictly on volume effects by using Eq . (4) 

(point 2). This adjustment is done independently of 
the stress levels. Because stress is not involved, 
elemental life is not affected, and the adjustment 
results in purely vertical displacement from the 
reference point. In the example shewn, the second 
element's volume is greater than the critical volume 
and contains more potential crack initiation aites. 
Therefore, the adjusted survivability ia lower. 


The element's life is then adjusted strictly on 
the basis of elemental stress relative to the critical 
stress by using Eq. (2> (point 3). This adjustment ia 
done independently of volume considerations. Conse- 
qusntly, elemental survivability is not affected, and 
the life adjustment results in purely horizontal dis- 
placement at constant survivability. In the example 
shown, the second element's stress is lower than the 
critical stress. Therefore, its life is greater and 
shifts to the right. 

A new material Weibull line ia now drawn through 
the adjusted life/survivebility point. Because the 
compared element's material is the same as the criti- 
cal element's material, the line is drawn parallel to 
the original line. However, if the elements had dis- 
similar materials, the new line would be drawn with 
the new Weibull slope. Thus, designs incorporating 
different materials can be compared with respect to 
survivability . 

By using Eq. (7), the compared eleinent's surviv- 
ability is adjusted once more to the critical ele- 
ment's assigned life of unity (point 4). Thin «■ j me 
the displacement is along the Weibull line and gives 
the survivability of an element relative to the criti- 
cal element, with both elements having the name life. 
Repeating this procedure for all the elements in the 
finite element model gives the survivability for any 
section of the design at any time. 

It should he mentioned that c burnt - Ve -I elemen- 
tal survivability will be somewhat in Tusnceo bv the 
finite element mesh, especially with components that 
experience fairly uniform 6tres?. For this case 
volumetric ratios will have more effect on elemental 
survivability. For components with nonuniform stress 
distributions, exponential stress irt icf will have 
more influence on elemental survivability. For tn^E 
case the methodology complements traditional finite 
element modeling philosophy, in that high-st resc - 
gradient areas should be modeled with smaller ele- 
ments. In this way sizing effects are minimized. In 
both cases an adaptive meshing techniu ie could be used 
to explicitly optimize the finite element model sizing 
and to implicitly minimize the mesh size effects on 
elemental survivability. From Eos. (4) and (8) it 
can be shown that component survivability remains 
unaffected by the individual element mesh. 

RESULTS AND DISCUSSION 

Parametric analytical Btudies were conducted to 
calculate speed effects on the survivability of a 
generic disk and to examine the probabilistic effects 
of material properties on disk survivability. The 
generic disk represents the first step in the investi- 
gation of aerospace propulsion turbine disks. It is 
61 cm (24 in.) in diameter and has \welve 1 . 02 -cm 
(0.40-in.) diameter bolt holes. 

Figure 2(a) shows the finite element model of the 
disk used in these studies. From the ax isymmetr ic 
conditions of th ' disk a IS degree Bectior of the disk 
is modeled with b2 grid points and 42 eight noded 
solid elements. Boundary conditions, defined in verms 
of the model's cylindrical coordinates, constta.r, dis- 
placement in the circumferential direction and a 1 ow 
displacement in the radial and axial directions. The 
disk stresses are a result of the centrifugal loads 
due to disk rotation. No thermal loads were included. 

MSC/NASTRAN was used fot the finite element code. 
The linear elastic analysis option, SOL 24, was used 
to calculate the stresses. Although this analysis is 
limited to stress levels below the yield point, it can 
be used for high cycle fatigue where stresses are 
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low«r and the occurrence of failura la laas determin- 
latic and more probab ilia tic. Howavar, bacauaa the 
statistical failura thaory ia baaad on atraaa lavala 
and dietribution, tha lifa and raliability methodology 
can alao ba adaptad for nonlinaar analyses. 

Tabla I liata tha assumed iaotropic matarial pro- 
part iaa uaad in tha analyaaa. Tha probabiliatic 
aapact of tha material strength ia given in tha L 
lifa of tha matarial. Thia L J0 lifa rapraaanta tha 
lifa in tha numbar of atraaa cyclaa in which 90 par- 
cant of tha taatad coupona aurviva at a conatant maxi- 
mum ahaar atraaa. By uaing tha valuaa of Tabla I, a 
finita alamant atraaa analyaia waa parformad on tha 
generic disk of Fig. 2(a) at a spaad of 12 800 rpm. 
Figura 2(b) ahowa tha raaulta of tha atraaa analy- 
aia. Tha maximum ahaar atraaa waa approximataly 
275 800 kN/m 2 (40 kei) at tha bolt hola. Tha failura 
probability analyaia for aach of tha modal' a alamanta 
waa then calculatad by uaing Eq. (7). Tha raaulta 
ahown in Fig. 2(c) indicata that tha region with tha 
lowaat raliability <i.e., the highaat probability of 
failure) occurred at, or near, tha maximum ahaar 
st rasa region at the bolt hola. Howavar, outaide thia 
region tha diak had relatively high ragiona of 
survivabil ity . 

Figura 3 graphically ahowa tha relation between 
lifa and aurvivability for tha matarial in Tabla I at 
an aaaumad Weibull slope. For a Waibull elope of 
3.57, 99 percent of the samples survived 6000 cycles, 
90 percent survived 10 000 cycles, and 10 percent sur- 
vived 20 000 cycles. A higher value of Weibull slope 
would demonstrate leas scatter, with failures occur- 
ring over a smaller number of stress cycles. A lower 
value would show failures over a broader range of 
Uvea and therefore might be indicative of more 
uncertainty in tha material lives. 

Figure 4 shows the effect of a disk's rotational 
speed on its probability of achieving a life of 10 000 
sties* cycles. The speed effect was plotted for an 
assumed Weibull slope of 3.57 and a strass-lifa axpo- 
nent of 9. The disk first experienced a maximum shear 
stress of 275 800 kN/m : (40 ksi) at a spaed of 
12 800 rpm. This disk had a probability of survival 
of about 96 percent at this speed, which waa higher 
than the coupon survival rata of 90 percent at the 
same stress and number of cycles. 

The higher probability of survival waa due to the 
volumetric effects. Tha 275 800-kN/m 2 (40-ksi) stress 
region in tha disk was much smaller than tha gaga vol- 
ume of tha coupon. The location of the maximum stress 
in tha model occurred at the elements adjacent to the 
bolt hole, at the three o'clock position. Alao, the 
stresses in tha remainder of tha diak elements were 
much lower than 275 800 kN/m : (40 ksi) and had individ- 
ually high survivability. Therefore, they had hardly 
any effect on degrading the disk's survivability. As 
a result, the disk showed a probability of survival 
greater that the coupon L <c life. 

The individual calculated elemental survivabili- 
ties can also be examined for design opt imi tat ion . 
Regions with low-survivability elements can ba rede- 
signed to increase reliability. Similarly, regions 
containing elements with unnecessarily high surviva- 
bility can be optimized through redesign to reduce 
weight . 

Tigure 5 ahowa the effect of apeed on tha cumula- 
tive probability of failure (i.e., 1 minus tha aurviv- 
ability) with varying Weibull slopes. Thia figure 
demonstrates how the distribution of failure ae a 
function of Weibull slope affects the probability 
of survival. All three curves intersect at 
12 800 rpm, with a 4 percent failure (96 percent 
survival), because the life reference street of 

275 800 kN/m z (40 kai) firat happens at 12 800 rpm. 
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Above thie apeed the rate of failure was definitely 
higher with higher Weibull alope. The 50 percent 
failure probability occurred at 13 100 rpm for a 
Weibull alope of 2, at 13 300 rpn for a Weibull slope 
of 3.57, and at 14 000 rpe for a Weibull slope of 5. 

Below 12 800 rpm the value of the Weibull alope had 
little effect on survival probability. 

Often, material data will be incomplete, if 
available at all, and designers must use engineering 
judgment in estimating material uncertainties. Where 
Weibull parameters are unknown for a rotating diak, it 
may be permissible to assume the Weibull alope and tha 
atrees-life exponent, in order to determine the diak 
speed at which the probability of survival will be 
high. From Tig. 5, the esv.imate of the Weibull alope 
should be high for conservative design by anticipating 
a higher failure rate than would probably occur. 

Figure 6 ahowa the effect of speed on the cumula- 
tive probability of failure with varying stress-life 
exponents. There is not the degree of divergence of 
the three curves above 12 800 rpm as there waa with 
Weibull alope variations. The same general trends 
hold, higher failure rates with higher atreaa-life 
exponent, but not to the same extent as in Fig. 5. 

At high survivability values, the value of the stress- 
life exponent had little effect at any given apeed. 

CONCLUDING REMARKS 

Probabiliatic material properties expressed in 
terms of Weibull parameters were coupled with the 
stress field determined from* MSC/NASTRAN finite 
element analyaia to determine fatigue life baaed on 
crack initiation. Thia methodology can be applied to 
predict the probability of survival of the complete 
structural component aa well as to identify critical 
failure regions of the component. A unique advantage 
of thia approach ia that only coupon fatigue testing 
is needed to establish the material fatigue parameters 
necessary for full-size componert life and survivabil- 
ity analysis. Thus, thia technique can be used in the 
c?rly design stages to optimize life-based designs, 
thereby reducing the amount of full-size component 
testing required to validate these designs. 

i 

SUMMARY OF RESULTS 

A method for calculating a component’s design 
survivability that incorporate* finite element analy- 
sis and probabiliatic material properties waa devel- 
oped. The method evaluates design parameters through 
direct comparison of component survivability expressed 
in terms of Weibull parameters. The method considers 
the component ' a total stress state in the survivabil- 
ity calculations. Critical regions with respect to 
aurvivability can be identified for optimization 
purposes: Reliability can be improved at low- 

aurvivability regiona; and weight can be reduced in 
high-survivability regions. Tne analyaia waa applied 
to a rotating disk with mounting bolt holes. The 
following results were obtained: 

(1) The highest probability of failure for the 
diak occurred at, or near, the maximum shear stress 
region at the bolt hole. 

(2) Diatribution of failure aa a function of 
Weibull alope affected the probability of survival. 

For speads that induced stresses abov# the L lifa 
reference stress, highar Waibull alope predicted lower 
survival probability. However, tha value of the 
Weibull slope had little effect at any given apeed at 
high probabilities of survival. 
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(3) The etrees-life exponent effected the 
reliability prediction*. Tor speed* that induced 
• tr ®**** *bove the L lQ life reference stress, 
higher stress-life exponent* predicted lower survival 
probability. However, the value of the strees-life 
exponent had little effect at any given speed at high 
probabilities of survival. 

(4) Where Weibull parameters are unknown for a 
rotating disk, it nay be permissible to assume Weibull 
parameters as well as the stress-life exponent in 
order to determine the disk speed at which the 
probability of survival will be highest. 
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TABLE I. - ASSUMED ISOTROPI C MATERIAL PROPERTIES 
Elastic modulus, kN/m : (pei; . . 1.10X10 11 (16X10 6 ) 


Poisson's ratio C.33 

D«n»ity , kg/m 3 (lb/in. 3 | 4.61*10' 5 (0.16) 

G«g« volum., m 3 (in. 3 ) 5.74X10' 6 (0.35) 

* l0 number of cvcles 10 000 

Shear stress at L, c , kN/m : (ksi) . . 275 800 (40) 



Figure 1 —Element survivability methodology. 
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(a) Disk finite element model. 
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(b) Shear stress contours at detail A. 


Life, number of stress cycles 

Figure 3.— Assumed material survivability at 
275 800 kN/m? (40 ksi). Weibull slope. 3 57 
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(c) Probability of survival contours at detail A. 

Figure 2.— Disk finite element model and analysis 
results. Nominal disk configuration: disk diameter, 
0 061 -m (24-in ); disk speed, 12 800 rpm; twelve 
1 .02-cm (0.40-in.) diameter bolt holes located on 
0.030-m (12-in.) diameter. 




Figure 4. —Effect of speed on surviv 
ability at 10 000 cycles Stress-life 
exponent, 9; Weibull slope, 3 57 







Report Documentation Page 


1. Report No. 

NASA TM -104400 
4. TWaand SuMMa 


2. Qovammao* Aeoaaaion No 


Incorporating Finite Element Analysis Into Component Life 
and Reliability 


3. Radptenfs Catalog No. 


5. Report Dele 


$. Performing Organization Code 


7. Authors) 

Richard August and Erwin V. Zaretsky 


9. Performing Organization Nairn and Addraaa 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

12. Sponaoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 - 0001 


«. Performing Organization Report No 
E -6220 

10. Work Unit No. 

505-63-5B 

11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Memorandum 


|14 Sponsoring Agency Coda 


1g Cli rTVrrrmtnry Nntm ■ - I 

Prepared for the Ninth Biennial Conference on Reliabiliiy. Siress Analysis, and Failure Prevention snon-ored by the 
Mechanical Engineers, Miami, Florida, September 22-25, 1991. Richard AugtS, Sverdrup 

IaS^S^NAsT 25 STf V Ae ™f“ P " kwa >'' Bt °° k P " k - 0h ‘» «'« (work Jdcd by 

Z " 0,Sk> ' NASA *«“**> C,™« Responsible pemon, RlchaM 

1ft Abelrerd " 1 

A method for calculating a component's design survivability by incorporating finite element analysis and n-obabi- 
'C P™*™ 5 »“ n* "Kthod evaluates design paranfetms through dnm^»ri£f„f 

eomponenlsurvtvabillty expressed in terms of Weibull parameters. The analysis was applied lo a tmaling disk wilh 
motmimg hoi, hote. Hk lughcsi pfobabiiiiy of failure occuned „, or near, iite maximl shea, sircss “fion afZ 
Distribution of failure as a function of Weibull slope affects the probability of survival. Whcfe Weibull 
D meters are unknown for a rotating disk, it may be permissible to assume Weibull parameters, as well as the 
stress-life exponent, in order to determine the disk speed where the probability of survival is highest. 
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